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ABSTRACT 

We introduce a likelihood analysis of multi-epoch stellar line-of-sight velocities to constrain the 
binary fractions and binary period distributions of dwarf spheroidal galaxies. This method is applied 
to multi-epoch data from the Magellan/MMFS survey of the Carina, Fornax, Sculptor and Sextans 
dSph galaxies, after applying a model for the measurement errors that accounts for binary orbital 
motion. We find that the Fornax, Sculptor, and Sextans dSphs are consistent with having binary 
populations similar to that of Milky Way field binaries to within 68% confidence limits, whereas the 
Carina dSph is remarkably deficient in binaries with periods less than ~10 years. If Carina is assumed 
to have a period distribution identical to that of the Milky Way field, its best-fit binary fraction is 
0.141q q5, and is constrained to be less than 0.5 at the 90% confidence level; thus it is unlikely to 
host a binary population identical to that of the Milky Way field. By contrast, the best- fit binary 
fraction of the combined sample of all four galaxies is 0.46^0 Jg, consistent with that of Milky Way field 
binaries. More generally, we infer probability distributions in binary fraction, mean orbital period, 
and dispersion of periods for each galaxy in the sample. Looking ahead to future surveys, we show 
that the allowed parameter space of binary fraction and period distribution parameters in dSphs will 
be narrowed significantly by a large multi-epoch survey. However, there is a degeneracy between 
the parameters that is unlikely to be broken unless the measurement error is of order ~0.1 km/s or 
smaller, presently attainable only by a high-resolution spectrograph. 
Subject headings: binary stars, theory — galaxies: kinematics and dynamics 



1. INTRODUCTION 

Recent progress in hydrodynamical simulations has al- 
lowed, for the first time, numerical simulations of star 
formation that include all the re l evant physics dow n 
to solar system scales (iBatd 120091 : lOffner et all 120091 ). 
Star formation is now understood to occur due to the 
gravitational collapse and fragmentation of a turbulent 
molecular cloud, gen erally in the presenc e of radiative 
feedb ack (|Batell201^ and magnetic fields (jPrice fc Bate! 
I2010r ). Simulations of star formation that incorporate 
these physical effects now offer detailed predictions about 
the statistical properties of stellar systems that must be 
tested against observational data to arrive at a complete 
theory of star formation. A major component of such 
predictions involve the statistical properties of binary 
star systems. While approximately half of solar-type 
field stars in the Milky Way are known to have binary 
companions, the fraction of young pre-main sequence 
stars in binary or higher-order star syst ems is much 
higher and possibly greater than 90% (|Leinert et al.l 
fimlKoider fc Leinertl(l99llPatience et al.ll2002r ). This 
binary fraction is later reduced as many wide binaries are 
subsequently disrupted by cluster dynamics or orbital de- 
cay. Thus, binary systems are likely the dominant mode 
of star formation, and the ability to successfully predict 
their statistical properties will be an essential component 
of a complete theory of star formation. 

At present, a complete census of the statistical prop- 
erties of binary star systems exists only in the solar 
neighborhood, at dist ance s out to » 30 pc from th e 
Sun (.Raghavan et al..,20Tol iDuquennov fc Mavoi]ll991[ ). 
While a number of visual binary studies have been 



made in ne arby open clusters (jBrandner fc Koehlei]ll998l : 
iPatience^t al. 2002), these studies are limited to wide bi- 
naries with separations > 10 AU. Moreover, such wide 
binaries are likely to have been "dynamically processed" 
by stellar encounters within the cluster, and thus may 
not represent the primordial binary populations present 
shortly after star formation has occurred. Likewise, spec- 
troscopic and p hotometric s tudies have been made in 
globular clusters (lYan fc Cohen 19 96: Milone et al. 20121 
iSollima et al.l 120071: iSollima et al.ll2012l ). but in many of 
these systems even close binaries have been dynamically 
processed in the dense core of the cluster over cosmo- 
logical timescales. To test theories of star formation, it 
is important to study binary populations in diffuse clus- 
ters or field populations over a range of separations for 
which the binary systems can be considered primordial, 
and compare their statistical properties to that of simu- 
lations. 

Dwarf spheroidal (dSph) galaxies are excellent objects 
for this purpose. Because they are more diffuse than 
large globular clusters, stellar encounters are sufficiently 
rare in dSphs that binary stars with separations less than 
10 AU should be largely unaltered since their formation. 
Dwarf spheroidals represent a variety of star forming en- 
vironments, ages, and metallicities, against which theo- 
ries of star formation can be tested. Furthermore, as sim- 
ulations advance and our understanding of star formation 
is refined, dSphs may ultimately become laboratories for 
a new brand of "stellar archaeology" : using the statisti- 
cal properties of their binary star populations to discern 
the initial conditions under which the galaxy originally 
formed. The extent to which binary properties are sensi- 
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tive to initial conditions during star formation is still an 
open question. While differences between bin ary popula- 
tions of different clusters have been observed (|King et alj 
l2012f l. these differences may be largely due to subsequent 
dynamical processing with the cluster (jMarks fc Kroupal 
12012 . ,Kroupa,1995 ). If binary star properties do indeed 
have a significant dependence on initial conditions (e.g. 
cloud temperature, density, magnetic field, metallicity), 
they may eventually become a tool for understanding 
star formation histories. 

Another reason for constraining binary properties in 
dwarf galaxies is to understand the effect of binary or- 
bital motion on the stellar velocity distribution, and 
hence the estimated mass distribution, o f the galaxies 
(jOlszewski et al.lll996t iMinor et al']|2010D . The spatial 
mass distribution of dark matter within a galaxy, includ- 
ing the existence and size of a central core, is determined 
by properties of the dark matter particle itself. How- 
ever, because there e xists a degeneracy between mass and 
velocity anisotropy (jWolf et al.ll2010( ). it is not possible 
to infer the dark matter density profile in a straight- 
forward way. One way to overcome this degeneracy is 
to make use of higher moments in the velocity distri- 
bution to distinguish between radial and t angential or- 
bits, and hence infer the velocity anisotropy (iLokas et ahl 
[200a lLokasl[2009t IRichardson fc Fairbairnll2012D . How- 
ever because binary motion can also contribute signif- 
icantly to higher moments in the velocity distribution, 
it is important to correct for binaries when determining 
the velocity anisotropy. Therefore, quantifying the effect 
of binary motion on a galaxy's velocity distribution can 
lead to a better understanding of the distribution of dark 
matter in dwarf galaxies. 

Finding detailed constraints on the binary populations 
in these galaxies is a difficult task. In the most sensitive 
spectroscopic surveys of dwarf spheroidals to date, typi- 
cal line-of-sight velocity measurement errors (« 1-3 km/s 
for red giant stars) are large enough to render binary mo- 
tion in systems with periods longer than 10 years unob- 
servable. Since binaries with such long periods cannot be 
directly observed in these galaxies, any statistical infor- 
mation about long-period binaries can only be inferred by 
extrapolation. Furthermore, among the red giant popu- 
lation, information about the shortest-period binaries is 
lost because close binary systems are d estroyed when th e 
red giant star overffows its Roche lobe (|Paczyhskilll971l ). 
Another complication is that any component of a star's 
velocity due to binary motion must be separated from the 
measurement error, and therefore deriving accurate and 
reliable measurement errors is essential before obtaining 
binary constraints. 

In this paper, we show that all of these difficulties can 
be surmounted by a likelihood approach, and demon- 
strate the method using multi-epoch data fro m the Mag- 
ellan/MMFS sample of I Walker et"all (|2009aD in the Ca- 
rina, Fornax, Sculptor, and Sextans dwarf spheroidal 
galaxies. In section [2] we derive a multi-epoch likelihood 
for binary and non-binary stars, followed by a discussion 
of binary model uncertainties in section [3l In section 
m we show how this likelihood can incorporate an er- 
ror model to derive accurate measurement errors in the 
Magellan/MMFS sample, which we apply to the data in 
section [5j To ensure the error model has properly repro- 
duced the measurement errors, we then test the model 



on simulated data in section [51 Next, the binary fraction 
in each galaxy is inferred in section [3 under the assump- 
tion of a Milky Way-like period distribution, and our 
main results are plotted in Figure IH wherein probabil- 
ity distributions in binary fraction are plotted for each 
galaxy. The inferred binary fraction of the combined 
sample of all four galaxies is plotted in Figure [5l and the 
best-fit binary fractions in each galaxy are listed in Ta- 
ble [2] In section |8] we find more general constraints on 
the binary fraction, mean period, and spread of periods 
in each galaxy. In section [9] we apply our methodology 
to simulated data to investigate the prospects for finding 
more precise binary constraints in a larger multi-epoch 
dataset. Finally, we discuss implications of our results in 
section [10] and conclude in section [TTJ 

2. MULTI-EPOCH LIKELIHOOD 

In order to constrain binary properties from radial ve- 
locity measurements, we must find the likelihood of our 
binary model parameters for a given set of measured ve- 
locities. We take as a model parameter the fraction of 
stars in binary systems B. Then the likelihood will take 
the form C = Lnb + BCh, where and Cb represent the 
likelihood for non-binary and binary stars, respectively. 
We shall further model the binary population by a set of 
parameters V that characterize the distributions of bi- 
nary properties. In general, these binary properties may 
include the periods, mass ratios and orbital eccentrici- 
ties. The distribution of these properties and our choice 
of model parameters will be discussed in detail in section 

In principle, even in a sample of stars with only one 
velocity measurement each, binary properties could be 
inferred from the line-of-sight velocity distribution of the 
sample, wherein the presence of binaries is manifest as a 
non-Gaussian tail representing short-period binaries with 
large orbital velocities. However, inferring binary prop- 
erties from the velocity distribution is fraught with dif- 
ficulties. Because the classical dwarf spheroidal galaxies 
have velocity dispersions greater than 7 km/s, only stars 
with radial velocities that differ from the galaxy's sys- 
temic velocity by at least > 14 km/s will be evident as 
binary stars, and only a small fraction of binaries have 
velocity amplitudes this large. For this reason, extremely 
large samples would be required to have any hope of 
constraining binary properties from single-epoch veloc- 
ity measurements alone. Contamination from Milky Way 
stars must also be accounted for, as such stars may be 
mistaken for binaries. Finally, non-Gaussianities may be 
present in the galaxy's velocity distribution due to ve- 
locity anisotropy and must also be distinguished from 
binary motion. 

Because the present work is primarily concerned with 
finding binary constraints, we can avoid the difficulties 
inherent in using the galaxy's velocity distribution by 
focusing on the stars with repeat measurements. Sup- 
pose a star has a set of n velocity measurements {vi\ = 
{ui, • • • ,«„} and errors {ui} taken at the corresponding 
dates, {ti\. For stars with at least two measurements 
{n > 2), we shall derive a likelihood in velocity differences 
Awi, A7J2, and so on, where we define Aw^ = Vi+i — vi. 
For readability, when denoting probability distributions 
we will suppress the brackets denoting sets of measure- 
ments (e.g., P{{vi\) P{vi)). 
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Our desired likelihood is derived as follows. We first 
write down the probability, for single stars, of drawing 
a set of velocity measurements {vi} given that it has in- 
trinsic velocity Vtme- Assuming the measurement errors 
are Gaussian, this likelihood can be written as follows: 



= 7V(t;.,,CTi)- 



(1) 

Here, (v) and am are the variance-weighted average ve- 
locity and measurement error. 



9 ' 



(2) 



(3) 



and a tedious algebraic calculation shows that the factor 
Af{vi,ai) is given by 



J\f{vi,ai) 



X exp ■ 



1 

4 ^ 



« 3 t J 



(4) 



The last term in the denominator of the exponent is ex- 
plicitly zero when n — 2. Note that this factor depends 
only on velocity differences Vi — Vj ; it has no dependence 
on the star's intrinsic velocity vtme- 

To derive a likelihood in velocity differences Aw, 
we change variables in eq. [1] from {vi^V2,v^, ...) to 
((u),Awi,Au2, ...) and integrate over (w). Then the 
Gaussian factor in eq. [1] integrates to 1, and we are left 
with 



ulated stars and binning in the velocity differences Af^ 
(t he details of the M onte Carlo simulation are discussed 
in lMinor et al.|[201ClD . The likelihood depends not only 
on the velocities and measurement errors, but also the 
time intervals between epochs and the magnitude M of 
the star, from which we determine the mass and radius of 
the primary star via a stellar population synthesis model. 
The likelihood will also depend on the set of model pa- 
rameters V used to characterize the binary population, 
which we discuss further in section [3l In terms of all the 
relevant model parameters, the likelihood for each star is 
therefore 



C{Av,\(y,,U,M-B,V) 

= {l- B)£„6(Aw,|a,) + BCb{^v,\a^,U,M■V). (6) 

We can put this in a somewhat simpler form by noting 
that the likelihood for single stars does not depend on 
any model parameters, and thus can be divided out of 
eq. m We then have 



C{^v,\a,,U, M; B, V) {I - B) + BJ{V), (7) 



where 



Cnb{Avi\ai) 



(8) 



For non-binary stars or binaries with small velocity 
amplitudes compared to the measurement errors, the J- 
factor in eq. [7] will be of order 1 or less, since the bi- 
nary likelihood will typically be slightly smaller than the 
non-binary likelihood. However, if a binary star exhibits 
large velocity changes compared to the measurement er- 
rors, the J- factor will be larger than 1, possibly by many 
orders of magnitude, since the non-binary likelihood Cnb 
in the denominator will be very small. 

For each star, the factor J('P) can be calculated over a 
grid of values in the binary parameters 7^, after which the 
J-factors can then be found by interpolation for any val- 
ues of the binary parameters V. After calculating the J- 
factors, the binary fraction B and parameters V can then 
be inferred either by a maximum-likelihood or Bayesian 
analysis. 



Cnb{^Vi\(7i 



-1 



^ ^ X ^ ~ 



(Aw, 



(5) 



This is our desired multi-epoch likelihood for non- 
binary stars; it is essentially the same factor Af{vi,ai) 
above, but expressed in the new variables Avi. As an 
aside, we note that it is possible to derive a likelihood in 
all the velociti es Vj using this formah sm; this is the ap- 
proach taken in [Martinez et al] (j201lD . Since we are only 
interested in binary constraints rather than the galaxy's 
intrinsic velocity distribution, only a likelihood in veloc- 
ity differences is considered in the present work. 

The corresponding likelihood for binary stars cannot 
be found analytically, so it must be calculated by run- 
ning a Monte Carlo simulation for a large number of sim- 



3. BINARY POPULATION MODEL UNCERTAINTIES 

Here we address the question of which binary model 
parameters to vary and how to deal with uncertainties 
in these parameters. Besides the binary fraction, a pop- 
ulation of binary stars can be described by distributions 
in three parameters: the mass ratio q, eccentricity e, 
and orbital period P. In the absence of a large number 
of epochs, eccentricities are difficult to constrain because 
very eccentric binaries spend a relatively small amount of 
time near their perihelion where the observed velocities 
are large. We therefore fix the distr ibution of ec centric- 
ities and assume the form given in iMinor et al.l ([2010), 
which is similar to that observed in solar neighborhood 
field binaries. 

Along similar lines, velocity measurements at several 
epochs are usually needed to determine the mass ratio 
of a binary independently of its orbital period. Theo- 
retical considerations suggest, however, that the period 
distributions of different binary populations may differ 
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drastically (|Brandner fc KoehleHll998l ). while the distri- 
bution of mass ratios may have a more nearly universal 
form. This is certainly true for long-period binaries, for 
which the mass ratio follows the Salpeter initial mass 
function for q > 0.5 (assuming the primary mass to lie 
in a restricted ran ge, which is the case for s tars older 
than a few Gyr; cf. iDuauennov fc MaYoHll99ll ). The ob- 
served distribution of mass ratios for sh ort-period bina- 
ries (P < 1000 days] is clo ser to uniform (jGoldberg et al.l 
120031 : iMazeh et al.lll99^ . and at present it is unclear 
whether this form is universal in primordial binary pop- 
ulations. We therefore fix the mass ratio distribution 
and assume it to follow a form similar to t hat observed 
in the solar neighborhood, as described in iMinor et al.l 
(|2010f ). with a uniform distribution for short-period bi- 
naries. Note that we are allowing for the mass ratio and 
eccentricity to vary from star to star — it is just the form 
of the distribution from which these parameters are de- 
rived that is fixed. In principle, we could also vary the 
functional form of the mass ratio distribution, but this 
is computationally expensive. The main reason is that 
the binary likelihoods will have to be computed on a grid 
that includes the parameters used to describe the mass 
ratio distributions. In addition, without a large number 
of epochs to separate each star's binary period from its 
velocity amplitude, these parameters will be highly de- 
generate with other binary parameters (specifically, the 
binary fraction and period distribution parameters). 

For simplicity we assume the period distribu- 
tion of dwarf spheroidal galaxies to have a log- 
nor mal form, in analogy to Milky Way field bina- 
ries (IDuauenno v fc Mavoi^l l99lUFischer fc Marcvlll992l: 
iMavor et aLr fT992: Raghav ajTet al.l |2010[) . Initiallv we 
shall fix the mean period fiiogp and spread of periods 
CTiogp to the values observed in solar neighborhood field 
binaries (where fJ.iogP — 2.24, txiogp — 2.3), while in later 
sections dHE]) these parameters will be allowed to vary. In 
general, we therefore have three binary parameters that 
will be constrained: the binary fraction B, the mean log- 
period /iiogP, and log-spread of periods ciogp. It must 
be borne in mind that for the dwarf galaxy samples con- 
sidered in this paper, binary motion in stars with periods 
> 10 years will be effectively hidden by the measurement 
error (the correspondence between observable binary pe- 
riods and measurement error will be explored in detail 
in section [9. 3p . Thus we can only directly probe the pe- 
riod distribution for periods shorter than ~10 years, so 
the shape of the period distribution at longer periods 
must be inferred by extrapolation. This extrapolation 
will hold only to the extent that the period distribution 
takes a symmetric log-normal form, which we take as a 
working hypothesis in this paper. 

4. MEASUREMENT ERROR MODEL 

Since any velocity variations beyond what is indicated 
by the measurement errors may be interpreted as binary 
motion, to constrain binary properties it is essential to 
estimate the meas urement e rrors a s well as possible. In 
the data from Walker et al.l (|2009a[ ). measurement errors 
are determined from repeat velocity measurements by 
applying an error model, which will be reviewed and ex- 
tended in this section. It is important to note, however, 
that if measurement errors are estimated by means other 
than velocity variability (e.g. from signal-to-noise prop- 



erties), an error model such as the one considered here 
need not be applied before obtaining bin a ry con straints. 

In the error model from IWalker et all (|2009aD . Gaus- 
sian measurement errors a are estimated using the model 



'ccn 



<7CCF 



(9) 



(10) 



where (Jq is a baseline error, accF is the spectrum cross- 
correlation error, and Ri is called the Tonry-Dav is R- 
value for the measurement (jTonrv fc Davis! 1l979[ ). de- 
fined as the height of the maximum peak in the spectrum 
cross-correlation function divided by the average peak 
height. Th is error model i s an e xtension of the error 
model from iTonrv fc Davis! ()1979() . which takes x = I. 
The MIKE spectrograph has two channels, "red" and 
"blue" , which operate over different wavelength ranges, 
so we must have two sets of error model parameters rep- 
resenting each channel. Our error model then has four 
parameters given by the set £ — {ored, ctUue, Xred and 
Xbiue}- The baseline error also depends on the channel 
color, and was determined to be cro,red = 0.6 km/s and 
ccbitie = 0.26 km/s. 

The i mplementation of th e error model described 
above in IWalker et al.l (j2007D has a few disadvantages. 
First, their method of constraining the error model pa- 
rameters in makes certain assumptions about the intrin- 
sic velocity of each star before deriving measurement er- 
rors. Since the star's intrinsic velocity cannot be well 
determined from repeat measurements without knowing 
the measurement errors in the first place, this may lead 
to bias in the error model. 

In contrast, since our multi-epoch likelihood is ex- 
pressed solely in terms of velocity differences, we can nat- 
urally incorporate the error model with no assumptions 
about intrinsic velocity. If the star's measurements all 
have time intervals of less than a week between them, bi- 
nary motion can be neglected and we can incorporate the 
error model into our likelihood by expressing the mea- 
surement error in terms of the error model parameters. 
Combining eqs. [Hj and llOi we can write the error in the 
ith measurement of the star as 



(1 + R^) 



2x^ 



(11) 



where Ci is the channel color of the ith measurement 
(either "red" or "blue"). This expression is then inserted 
into the likelihood from eq. [5j which we now write as 



C{^V^\R,,C^■,£) = 



nr=iN/2^ 



(12) 



Mor e importantly, the error model in IWalker et al.l 
(|2007f ) does not take binary motion into account, which 
may bias the errors if the time between measurements 
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is longer than a few weeks. We shall address this by in- 
corporating binarity into the likelihood of each star if at 
least one of the time intervals between measurements is 
long enough for binary velocity variation to be observed. 
In the process, we will constrain the binary fraction B 
and the error model parameters simultaneously. 

At this stage, the reader may be wondering why our 
binary model is being applied twice: first to find the 
measurement errors, then afterward to constrain the bi- 
nary parameters B and V. Why not constrain all the 
error model and binary parameters simultaneously? The 
answer is that this would be extremely computationally 
expensive, because the binary likelihoods must be recal- 
culated every time the error model parameters are varied. 
Fortunately, if we are (for the moment) only interested in 
deriving the best possible measurement errors, it is suffi- 
cient to assume a Milky Way-like period distribution and 
constrain only the binary fraction B in order to distin- 
guish binary velocity variation from measurement error. 
Furthermore, as will become evident shortly, the binary 
likelihoods become computationally expensive to calcu- 
late while varying the error parameters if a star contains 
more than one "long" time interval (> 10 days) between 
measurements. Thus, any extra measurements occurring 
at later time intervals longer than 10 days will be dis- 
carded for the purpose of determining the error model 
parameters, although they will be used later to help de- 
termine binary constraints. 

In summary, while it may appear cumbersome to ap- 
ply the binary model twice, the benefit of finding more 
detailed binary constraints after determining the mea- 
surement errors is that ultimately the binary parameters 
V will be constrained in addition to the binary fraction 
B, and all the available measurements will be used to 
make this determination. 

To incorporate binary motion into our likelihood, we 
need to be able to quickly calculate the binary likeli- 
hood for a given set of assumed measurement errors, 
which we continually update as the error model param- 
eters are changed. To accomplish this, we first run a 
Monte Carlo simulation to generate the binary likelihood 
in log I Aw I with zero measurement error, which we denote 
by Lb (log I Aw I; At, M) where M is the absolute magni- 
tude of the star and Ai is the time interval between mea- 
surements. The likelihood is generated in log | Awj , rather 
than Au, because the binary likelihood becomes singu- 
lar as Aw — 0. The binary likelihood is generated over 
a table of absolute magnitude values M and time inter- 
vals At. In practice the long time intervals are typically 
close to a multiple of one year, so for our purposes it is 
sufficient to pick time intervals of 1 year, 2 years, and 
so on. For each star with a long time interval between 
measurements, we interpolate in magnitude and choose 
the time interval (in multiples of one year) closest to the 
actual interval between measurements. 

Now suppose that the star has n velocity measure- 
ments, and the long time interval is between measure- 
ment k and fc -I- 1. To find the binary likelihood in 
the presence of measurement error, we must find the 
convolution of the binary likelihood and the likelihood 
from eq. [T^] (representing measurement error) in the ve- 
locity difference Aw^. By switching to an integral over 
log I Aw^ I , this can be shown to equal 



/:b(Aw,|i?„c,,At,;,Af;f) 
1 

= -j Lb(log|AwlJ;Atfe,M) 

x{/:(---,Awfc-|Aw^|,---|i?„c,;£)-H 
£(•••, Awfc + lAw^l, • • • c,- £)]d{\og lAw^l) 

(13) 

where £(Awi, • • • , Aw„_i £) is the likelihood in 
eq. [121 This is our binary likelihood expressed in terms 
of the error model parameters, which must be calculated 
by performing the integral every time the error model 
parameters are varied. By now it should be evident why 
only one long time interval {l 10 days) between mea- 
surements is allowed: multiple long time intervals would 
require multiple integrations to account for binary vari- 
ability over each interval, and this would be prohibitively 
expensive to compute while varying the error parameters. 

Finally, for each star with a time interval longer than a 
week, we write the likelihood in terms of binary fraction 
B as 

/:(Aw,|i?,,c„At,,M;B,£) 

= (1 - B)£„b(Aw,|i?„ c; £) + SA(Aw,|E„ c„ At,, M; £) 

(14) 

We therefore shall constrain the four error model pa- 
rameters £ and binary fraction B simultaneously. We 
emphasize again that once the measurement errors have 
been determined by this procedure, we shall go back and 
find more precise binary constraints by including all the 
measurements. The maximum likelihood error model 
parameters obtained by this procedure will be used to 
determine the measurement errors for each star before 
proceeding with the full binary analysis. 

5. APPLICATION OF ERROR MODEL TO DATA 

We now apply the error model described in the pre- 
vious section to derive measurement errors in the Ca- 
rina, Fornax, Sculptor, and Sextans dSph galaxies. Be- 
fore plunging into the analysis, we note that there are 
potential sources of non-Gaussian velocity variability be- 
sides binary orbital motion. The most common of these 
are due to false peak selection in the spectrum cross- 
correlation function, which occurs at low Tonry- Davis 
j?-value s and low signal-to-noise ratios. In lWalker et al.l 
(|2009al) . measurements with i?- values less than 4 are dis- 
carded for this reason. We take the additional step of dis- 
carding measurements with signal-to-noise smaller than 
1.2, which can also exhibit non-Gaussian error. This is 
due to the fact that Gaussian measurement errors are not 
reproduced properly if very low signal-to-noise measure- 
ments are included (the criterion for determining this will 
become clear shortly), indicating the likely presence of 
non-Gaussian error among those measurements. On the 
other hand, the derived measurement errors are not sig- 
nificantly affected for signal-to-noise thresholds greater 
than 1, except for cuts well above 2 when the sample 
size is reduced significantly. Furthermore, we find that 
measurements using the red channel with i?-values less 
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Fig. 1. — Cumulative distribution of Av/cr2e for each galaxy in the Magellan/MMFS sample of lWalker et al.l H2009al ). Here, Av = \v2 —vi 
is the difference between two successive line-of-sight velocity measurements and (T2e is the measurement error in Av, given by (Tj^ = ct^ +(7| 
where cri, a2 are the individual errors on each measurement. We plot the distribution using measurement errors determined from our error 
model (red solid line). To check that the data can indeed be reproduced by our derived error model parameters, we plot a distribution for 
simulated data using the same number of stars, derived measurement errors, and e pochs as in th e real data (dashed line). For comparison, 
the same distributions are plotted using the published measurement errors from [Walker et aTl t2009a ) (dot-dashed line). The expected 
distribution in the limit of a large sample without binaries is also plotted for reference (solid dark line). 

than 7 correlate strongly with low signal-to- noise (;< 1), 
due to the red channel's lower spectral resolution com- 
pared to the blue channel. Therefore, as a further pre- 
caution these measurements are also discarded to avoid 
non-Gaussian errors. 

When generating the binary likelihoods, the star's 
magnitude is used to infer its mass and stellar radius 
und er the assumption that it lies on the red giant branch 
(see iMinor et al.ll201Cll for details) ; thus it is essential to 
exclude horizontal branch stars from each sample. In 
the Carina and Sculptor samples we discard stars with 
magnitudes My > 20.3 to ensure that horizontal branch 
stars are not included. The Fornax and Sculptor samples 
do not extend to sufhciently faint magnitudes to include 
horizontal branch stars, so no magnitude cuts are neces- 
sary in those samples. 

As outlined in section |4l our error model presently al- 
lows only one "long" time interval (> 10 days) between 
measurements — any further measurements occurring at 
long time intervals are discarded, although they will be 



included later when deriving detailed binary constraints 
in sections [7] and El After making the aforementioned 
signal-to-noise, i?-value, velocity and magnitude cuts in 
the Fornax sample, there remain 458 stars with repeat 
measurements, the largest of any galaxy in the survey. 
Of these, however, only 201 stars have time intervals of 1 
year or longer between any pair of measurements. Many 
of the latter subset have measurements at three epochs, 
with time intervals of approximately a week and one or 
two years between measurements. Since none of the stars 
have multiple time intervals greater than 10 days, no 
further measurements were discarded for the purpose of 
determining measurement errors in Fornax. 

After applying our error model to the data, the re- 
sult s for the Fornax dwarf galaxy are plotted in Fig- 
For all repeat measurements in the sample. 



1(a) 



we plot the cumulative distribution of Av/<T2e, where 
At; = I W2 — I'll is the difference between two successive 
line-of-sight velocity measurements and cr2e is the mea- 
surement error in Av, given by crfg = (tJ+(t| where ai, oi 
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Galaxy 


Median error (km/s) 


Median error, published 


Carina 


2.32 


2.6 


Fornax 


1.71 


1.1 


Sculptor 


2.10 


1.7 


Sextans 


2.64 


2.6 



TABLE 1 

Median measurement error comparison 

are the individual errors on each velocity measurement. 
If the velocity changes are only due to Gaussian measure- 
ment errors, in the limit of a large sample this distribu- 
tion should approach the complimentary error function 
erfc [Az;/-\/2cr2e] (solid dark line). 

From the figure it is evident that Gaussian measure- 
ment errors are well-reproduced in Fornax for At; < 
1.2CT2e, comprising ~80% of the stars in the sample, 
while binary orbital motion associated with large ve- 
locity changes becomes manifest above this threshold. 
The maximum likelihood binary fraction in Fornax is 
B « 0.5, while the maximum likelihood error parame- 
ters are £ — {ared = 36.9, auue = 55.2, Xred — 1-06 and 
Xbiue = 1.37}. To test this, we simulated a dataset with 
the same derived measurement errors, epochs, and mag- 
nitudes as in the real data, and assigned new velocity 
values for each measurement assuming a star population 
with binary fraction B = 0.5. The di stribu tion for a typ- 
ical realization is plotted in Figure 1(a) (dashed line). 



and is quite consistent with the data given the scatter 
in samples of this size. By comparison, we plot the dis- 
tribution using the p ublished measurement errors from 
iWalker et al.l ()2009al) . from which it is evident that the 
measurement errors are underestimated even for small 
velocity variations. This discrepancy may be not only 
due to binaries, but also the presence of low signal-to- 
noise measurements biasing the error model. We thus 
conclude that the published measurement errors for the 
Fornax dSph are significantly underestimated; in fact the 
median measurement error in the sample using the pub- 
lished errors is 1.1 km/s, while the median measurement 
error from our model is 1.71 km/s, differing by a factor of 
~ 55%. The median measurement errors in each galaxy 
are listed in Table [TJ 

In the case of the Carina and Sextans samples, contam- 
ination by foreground Milky Way stars is a major concern 
because of the relatively low surface brightness of these 
galaxies. In addition, many stars selected for repeat mea- 
surements have velocities more than 3cr away from the 
galaxy's systemic velocity, where a is the galaxy's veloc- 
ity dispersion. Since binary velocity variations greater 
than 50 km/s are extremely unlikely, we therefore discard 
stars with velocities > 50 km/s away from the galaxy's 
systemic velocity as probable nonmember stars. In the 
Carina sample, there are 257 stars with repeat measure- 
ments; after making the rough velocity cut described 
above, together with a magnitude cut at My = 20.3 to 
get rid of horizontal branch stars, only 107 stars remain 
in the sample. Similarly, the Sextans sample contains 
203 stars with repeat measurements, among which only 
134 remain after making velocity and magnitude cuts. 
By contrast, the Sculptor sample contains 198 stars with 
repeat measurements, of which 190 remain after velocity 
and magnitude cuts, indicating relatively little contami- 
nation by nonmember stars. 

Applying our error model to the remaining sample in 



the Carina dSph, we obtain an extraordinarily low best- 
fit binary fraction, B = 0.05, the value of which will 
be refined in section [T) Using our derived measurement 
errors, th e resu lting distributions in Av/a2e are plotted 
in Figure 1 (b) Because of the relatively small size of the 
sample, tne distribution can differ significantly from the 
expected large-sample limit, but the absence of binary 
velocity variation is immediately apparent in the figure 
(we will return to this in detail in section [T]) . To test 
our model, we again simulate a dataset with the same 
derived measurement errors, epochs, and magnitudes as 
the real data, and assign new velocity values for each 
measurement assuming our best-fit binary fraction for 
Carina. Two random realizations are plotted; although 
the distributions can vary considerably because of the 
small sample size, we see that the actual distribution 
using our measurement errors is indeed consistent with 
the model. 

Finally, we apply our error model to the Sculptor and 
Sextans d Sph s amp les an d plot distributions in Av/a2e 



in Figures 1(c) and 1(d) 



For each galaxy, we plot dis- 
tributions using our derived measurement errors, along 
with distributions using the published measurement er- 
rors. From the figure it is evident that for both galax- 
ies, the published data underestimates the measurement 
error for large velocity differences, while in Sextans the 
error is overestimated for small velocity changes. By con- 
trast, the Gaussian error is well-accounted for in our de- 
rived measurement errors, with binary velocity variation 
showing up for velocity changes greater than w 1.4 km/s. 

While the principle goal of our analysis in this sec- 
tion has been to derive measurement errors, a rough es- 
timate of the binary fraction in each galaxy has been 
derived simultaneously. However, the binary constraints 
obtained in this section are only approximately valid, for 
several reasons: first, to save time during the likelihood 
calculation, for each star with a long time interval (> 
10 days) between measurements, the long time interval 
was rounded to a multiple of 1 year — thus the binary 
likelihoods were generated for multiples of 1 year with- 
out the need for additional interpolation. For some stars 
the rounding error is as long as three months, so the bi- 
nary likelihood is only approximately accurate in those 
cases. More importantly, as already discussed, certain 
repeat measurements were discarded if the star already 
contained multi-epoch measurements with one long time 
interval. Some of those discarded measurements may in- 
dicate binary velocity variation and thus are important 
for finding binary constraints. Finally, by inferring the 
binary fraction, we have assumed thus far that the dis- 
tribution of orbital periods in each galaxy is identical to 
that of Milky Way field binaries, which may not be the 
case. Ultimately, it would be preferable to constrain not 
only the binary fraction, but also the period distribution 
parameters in each galaxy. The full binary analysis will 
be performed in sections [7] and [51 

6. TEST OF ERROR MODEL ON SIMULATED DATA 

To verify that we have reproduced the Gaussian mea- 
surement errors well, we shall apply our error model to a 
series of simulated datasets with a distribution of Gaus- 
sian measurement errors si milar to that obtained in the 
Magellan/MMFS survey of IWalker et~all (|2009aD . This 
is accomplished by drawing a random Tonry-Davis R- 
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Aw/(T2e, where Aw = |w2 — wi| is the difference between 
two successive line-of-sight velocity measurements and 
<^2e = fi + erf where cti, (T2 are the in dividu al e rrors 



Fig. 2. — Cumulative distribution of Av/a2e for simulated data sets with 2-epoch measurements, with a time interval of 1 year between 
measurements. Here, Av = \v2 — fi | is the difference between two successive line-of-sight velocity measurements and a2e is the measurement 
error in Ad, given by o-jg = cr^ -|- cr| where cri, cr2 are the individual errors on each velocity measurement. Shown are two random 
reali zations of dat a sets with 200 and 500 stars, respectively, where the distribution of measurement errors is similar to those in the data 
from rWalker et "al . (2009a) and the absolute magnitude limit Myn^ = 1. In each figure we plot a distribution for a sample with no 
binaries, using measurement errors derived by applying our error model to the simulated data (solid red line), and a distribution for the 
same dataset using the "true" measurement errors of the simulation (dashed line); for comparison, we also plot the theoretical expected 
distribution in the limit of a large sample (solid dark line). We also plot distributions for a sample with a binary fraction oi B = 0.5 
and a Milky Way-like period distribution, again using the derived measurement errors from the error model including the effect of binary 
orbital motion (dot-dashed line) and true errors (dotted line) from the simulation. In all cases the distributions useing the true errors are 
reproduced well by the model. 

value for each measurement from a Gaussian distribution 
whose mean and spread are similar to those observed in 
the real data, while the spectrograph channel used to ob- 
tain each measurement is chosen to be either red or blue 
with equal probability. The "true" error model param- 
eters are chosen to be identical to those obtained from 
the Fornax sample (which have been obtained in section 
El). First we consider samples with measurements at two 
epochs, separated a year apart, with sample sizes of 200 
and 500 stars. In each case, we simulate a sample with a 
binary fraction of zero and apply our error model without 
correcting for binaries. We also simulate a sample with a 
binary fraction B = 0.5 and apply our error model with 
the binary correction. In all simulations, an absolute 
magnitude limit My = 1 is assumed. 

The results are displayed in Figure [5J As in the pre- 
vious section, we plot the cumulative distribution of 



on each velocity measurement. In Figures 2(a) and 2(b) 
we consider two typical random realizations of a 200 
star sample, where each star has measurements at two 
epochs separated by a year, and plot distributions using 
the "true" measurement errors for each simulation, as 
well as the derived measurment errors by applying our 
error model to the simulation. In each realization, the 
distribution using the derived errors follows the distri- 
bution using the "true" errors for both the 5 = and 
B — Q.^ populations; the percent difference between the 
derived and true errors is smaller than 10% for at least 
K, 80% of the stars in the sample, and smaller than 30% 
for w 95% of the stars in the sample. The Gaussian mea- 
surement errors are thus reasonably well-reproduced by 
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Fig. 3. — Cumulative distribution of l^vja^e for simulated data 
sets of 200 stars with 4-epoch measurements, where the first time 
interval between measurements is 1 year, and the remaining time 
intervals is 1 day. Plotted curves are identical to those in Figure [2] 
Note that the measurement errors are better reproduced compared 
to a simila r sa mple with measurements at only 2 epochs (compare 
Figs. |2(a)] and [2(b)| , although the binary fraction is not any better 
constrained because only one time interval is long enough for binary 
variability to be observed. 



the error model for most stars in the sample. For velocity 
changes Au/fT2e < 1, the distributions follow the large- 
sample limit, indicating the distribution of velocities be- 
low this threshold is Gaussian and binary behavior is not 
distinguishable from measurement error for s uch sm all 
veloc ity changes. As can be seen from Figures 2(c) and 
2(d) the measurement errors are very well-reproduced in 
a 5 (JO- star sample regardless of the realization — in both 
cases shown, the percent difference between the derived 
and true errors is smaller than 10% for at least 97% of 
the stars in the sample. Since the distributions are cu- 
mulative, if a few measurement errors are reproduced 
incorrectly, the distribution using the derived errors can 
be offset from the true distribution; note however that 
in fig ! 
distri 



2(c) the distribution for B = 0.5 parallels the true 



ibution, indicating most of the measurement errors 
are indeed accurate. 

As a final sanity check, we verify that our error model 
works for a larger number of epochs. In Figure[3]we con- 
sider a sample of 200 stars with measurements at four 
epochs, where the first time interval between measure- 
ments is equal to 1 year, and the remaining time in- 
tervals is equal to 1 day. Since time intervals shorter 
than a week are too short for binary orbital motion to 
be evident, we expect the additional epochs to aid in 
constraining the measurement errors, although without 
helping to constrain the binary fraction directly. We find 
that applying our error model to the 4-epoch sample re- 
produces the measurement errors very well, although the 
sample with B = 0.5 is not as well reproduced as in the 
500-star sample because the binary fraction is not as well 
constrained. 

7. BINARY FRACTION IN CARINA, FORNAX, SCULPTOR, 
AND SEXTANS DSPH GALAXIES 

7.1. Sample selection criteria; contamination by Milky 
Way stars 



We now adopt a Bayesian approach to infer properties 
of the binary populations in each galaxy. For simplicity, 
in this section we will assume the distribution of binary 
periods in each galaxy is identical to that of Milky Way 
field binaries (with a mean period of ~ 180 years and log- 
sprea d of periods equal to 2.3; see iDuauennoy fc Mavoil 
119911 ). an assumption that wiU be relaxed in section |51 
Since we have no prior knowledge of the binary fraction 
in dwarf galaxies, we will assume a fiat prior in B. As 
in the previous section, we make a magnitude cut to ex- 
clude horizontal branch stars in each galaxy, and also dis- 
card measurements with low signal-to-noise to avoid non- 
Gaussian errors that may bias the inferred binary frac- 
tion. As we will show later in this section, we find that 
our results are essentially unchanged for signal-to-noise 
thresholds higher than 1.2, except that the constraints 
become weaker for higher thresholds as a greater frac- 
tion of stars are being removed from the sample. To be 
conservative, we will discard measurements with signal- 
to-noise smaller than 1.4 in the calculations that follow; 
later we will investigate how the results differ as this 
threshold is varied. 

Another potential bias in the inferred binary fraction 
is contamination by foreground Milky Way binary stars. 
The majority of foreground stars are K-dwarfs in the 
Milky Way disk, many of which can be expected to have 
binary companions. We can expect larger velocity varia- 
tions in these binaries compared to binaries in the back- 
ground dSph galaxy, for two reasons: first, since the ob- 
served nonmember star lies on the main sequence, the 
secondary star may not be significantly dimmer than the 
primary star. Thus the measured radial velocity may be 
that of the less massive secondary star, for which larger 
velocity amplitudes are expected. Second, since the ob- 
served member stars lie on the red giant branch, very 
close binaries with high velocity amp litudes will hav e 
been destroyed by Roche- lobe overfiow (jPaczviiskill 1 9 7 ll) , 
whereas this is not necessarily the case for the foreground 
Milky Way stars since the observed star lies on the main 
sequence. Thus if there is significant contamination by 
Milky Way binaries, we can expect that it will most likely 
bias the inferred binary fraction to a high value. 

We can cut down on contamination significantly by 
performing a rough velocity cut as outlined in section 
[5l In addition, since the metallicity distribution is 
well-determined in each galaxy in terms of the Mg I 
triplet (A ~ 5170A) absorption line strength, we can 
use metallicity information to help distin guish between 
member and non-member stars. Following IWalker et al.l 
(2()Q9SI) , we model the distribution of Mg-triplet pseudo- 
equivalent widths (or magnesium strengths) for both the 
member and nonmember stars as a Gaussian, so that the 
member star population has mean magnesium strength 
w and dispersion in magnesium strength aw ; likewise, the 
nonmember stars have similar parameters wmw, <^w,mw- 
The values of each of these parameters characterizing the 
metallicity di s tributio n in each sample are determined in 
IWalker et al.l ([2009b). For each individual star we cal- 
culate the average magnesium strength w by performing 
a weighted average over all measurements w^, and like- 
wise we calculate the average measurement error in mag- 
nesium strength, Cw Assuming Gaussian measurement 
errors in w, we then define the likelihood in magnesium 
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Fig. 4. — Posterior probability distributions in bina ry fraction for 
each of the galaxies in the Magellan/MMFS sample of lWalker et al.l 
l|2009a!). where the binary period distribution is assumed to be 
identical to that of Milky Way field stars. In the Fornax galaxy 
(solid curve) the most probable inferred binary fraction is consis- 
tent with Milky Way field binaries, whereas the Carina dwarf (dot- 
dashed curve) has a very low inferred binary fraction, indicating a 
dearth of short-period binaries compared to the other galaxies. 



strength for member stars as 



exp 



-{w - 



(15) 



The corresponding Hkehhood for Milky Way stars 
is identical, except with the Milky Way magnesium 
strength distribution parameters wmw, <Jw,mw- For the 
member stars, we use the same velocity likelihood as in 
equation [Bl while for the nonmember stars we use the 
nonbinary velocity likelihood in equation [5j Assuming 
a fraction F of stars in the sample are members, our 
likelihood can be written 



C{l:\Vi,Wi\ew,i,(Ti,ti,M\B,V) (16) 
= (1 - F)LMw{^Vi\ui)Lw,Mw{wi\ew.i) 
+FC{Av,\a^,U,M;B,r)C^{w,\e^.^) 

Since we are not modeling the binary population of 
the Milky Way stars, there is still the possibility that 
a binary nonmember star may be counted as a mem- 
ber if it exhibits a large velocity change. However, the 
stellar metallicities in the Sculptor, Sextans, and Carina 
galaxies are significantly lower than that of the nonmem- 
ber stars, so the metallicity likelihoods will significantly 
reduce the impact of these stars on the inferred binary 
fraction. The Fornax dSph has a metallicity distribution 
similar to that of the nonmember stars, so the method 
outlined above cannot resolve nonmember contamination 
in Fornax; fortunately however, because of its high sur- 
face brightness, the fraction of nonmember stars in the 
Fornax sample is expected to be very low. 

7.2. Binary fraction constraints 

Adopting th e likelihood in eg. [TBI we use a n ested sam- 
phng routine ()Skillind 120041: iFeroz et al.ll2009D to obtain 
marginal posterior probability distributions in the binary 



Galaxy 


Binary fraction 


Carina 


^ -. . -4-0 9R 


Fornax 


^■0^-0.17 


Sculptor 


0-66_o.i9 


Sextans 


O.66I0 24 


Combined 


n 46+°-" 

"■^"-0.10 



TABLE 2 

Best-fit binary fractions, assuming Milky Way period 
distribution 



fraction B and member fraction F in each galaxy, assum- 
ing a flat prior in each parameter. With the exception 
of Sextans, we find the inferred member fraction in each 
galaxy is greater than 0.9 to within 99% confidence lim- 
its, indicating our rough velocity cut (discussed in sec- 
tionO successfully removed most of the Milky Way stars. 
In Sextans, however, the member fraction is constrained 
to be 0.75Q 7g, indicating contamination may still be an 
issue, although metallicity information can still distin- 
guish between member and nonmember stars in most 
cases due to Sextans' low mean magnesium line strength 
{w — 0.36A, compared to wmw — 0.84A for the fore- 
ground Milky Way stars). Since the member fraction in 
each galaxy is well-constrained by the metallicities in any 
case, the result is not strongly dependent on our prior in 
the member fraction F. 

In Figure m we plot posterior probability distributions 
in the binary fraction B after marginalizing over the 
member fraction F. The most probable inferred binary 
fractions in each galaxy are listed in Table [5J where the 
errors are given by the 15.87% and 84.13% percentiles of 
the probability distributions (corresponding to 1-cr error 
bars in a Gaussian distribution). 

To evaluate whether the inferred binary fractions in 
Figure S] are consistent with that of Milky Way field bi- 
naries, we note that binaries with G-dwarf primaries in 
the solar neig hborhood a re found to have a binary frac- 
tion of « 0.5 ()Duquenn ov & Mayor 1991); however, bi- 
nary fraction is known to be a function of mass, with 
smaller-mass prin iarics havi ng smaller binary fractions 
(|Fischer fc Marcvl [1992. Ra ghavan et al.ll2010D . In the 
solar neighborhood, iRaghavan et al.l ()2010f ) finds that 
primary masses in the approximate range of 0.75-1Mq 
correspond to an average binary fraction of ^0.4, while 
1-1. 4Mq primaries correspond to an average binary frac- 
tion of ^0.5. The galaxies in this study have multiple star 
populations with widely varying ages, and thus a range 
of primary masses on the red giant branch. While pho- 
tometry shows the Sculptor and Sextans dSphs have pre- 
dominantly very old s t ar populations w ith ages > 10 Gyr 
(|de Boer et al.ll2012aL [Lee et al.ll2009[ ). the Fornax dSph 
show s a wide spread of st ellar ages from 2 Gyr up to 13 
Gyr (|de Boer et al.|[2012b[) , while Ca rina's dominant stel- 
lar p opulation has ages of 4-7 Gyr (|Hurlev-Keller et all 
119981 ). On the red giant branch, this corresponds to 
stellar masses ranging from 0.8Mq up to approximately 
1.2Mq. We therefore assume an inferred binary fraction 
in the range 0.4-0.6 to be roughly consistent with that of 
Milky Way field binaries. 

By this standard, Fornax, Sculptor, and Sextans all 
have binary fractions consistent with Milky Way field bi- 
naries to within 68% confidence limits. Carina, by con- 
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Fig. 5. — Posterior probability distribution in binary fraction 
using the combined sample of galaxies in the Magellan/MMFS 
dataset, where the binary period distribution is assumed to be iden- 
tical to that of Milky Way field stars. Plotted are the combined 
sample of all four galaxies (solid curve) and the combined sample of 
Fornax, Sculptor and Sextans (dotted curve) without including the 
Carina dSph which has an anomalously low binary fraction. The 
most probable binary fraction over the entire combined sample is 
0.46j[g jg, consistent with the Milky Way field binary fraction of 
^ 0.5 for solar-type stars. 

trast, has a most probable inferred binary fraction of 
0.14, and its inferred binary fraction is less than 0.29 
to within 68% confidence hmits, and less than 0.57 to 
within 95% confidence limits. It must be borne in mind 
that here we have assumed Carina's period distribution 
to be similar to that of the Milky Way. More gener- 
ally, the results imply that Carina is largely devoid of 
binaries with periods (< 10 years) that would produce 
observable velocity changes over the relevant timescale 
of a few years. As will be made explicit in section [51 this 
could also be explained if the mean period of Carina's 
binary population is longer than that of the Milky Way, 
or has a smaller spread of periods, rather than having 
a very small binary fraction. Whatever the reason. Ca- 
rina's lack of short-period binaries appears to be incon- 
sistent with having a binary population similar to that 
of Milky Way field binaries. We will discuss the possible 
implications of this result for star formation in section 

m 

Continuing under the assumption of Milky Way-like 
period distributions in each galaxy, we can ask, what is 
the binary fraction of the combined sample of all four 
galaxies? There are 621 stars in the combined sample, 
affording much better constraints on binary fraction com- 
pared to the sample in each galaxy separately. The in- 
ferred binary fraction of the combined sample is plotted 
in Figured] (solid curve), and its best-fit value is given 
in the last entry of Table [H The most probable inferred 
binary fraction is 0.46, and the 68% confidence interval 
lies in the range 0.37-0.59, which is consistent with the 
Milky Way field binary fraction « 0.5 for solar-type stars. 
Finally, since Carina has an anomalously low apparent 
binary fraction, we also plot the combined sample of For- 
nax, Sculptor, and Sextans without including the stars 
from Carina (dotted curve). The most inferred probable 
binary fraction without Carina is 0.55, with the 68% con- 
fidence interval lying in the range 0.44-0.70, somewhat 



high but still consistent with Milky Way field binaries. 
7.3. Low binary fraction in Carina dSph? 

The apparent shortage of short-period binaries in the 
Carina dSph cannot be easily explained by sources of bias 
in the inferred binary fraction. Systematic errors such 
as contamination by foreground Milky Way stars, non- 
Gaussian errors due to low signal-to-noise or other fac- 
tors, or the presence of variable stars all tend to increase 
the non-Gaussian velocity variability of the stars, and 
thus would inflate the inferred binary fraction, rather 
than reduce it. Indeed, such biases might be a factor in 
the somewhat high inferred binary fractions of Sextans 
and Sculptor. In particular, contamination by Milky 
Way binaries could be infiating the binary fraction of 
Sextans, which has a lower member fraction (w 0.75) in 
our selected sample compared to the other galaxies. 

We can, however, investigate whether the Carina re- 
sults are sensitive to the assumed signal-to-noise cut- 
off. In Figure [6] we plot posteriors in the binary frac- 
tion of Carina for signal-to-noise cutoffs ranging from 1 
to 1.6. We find that increasing the cutoff from 1.2 does 
not change the inferred binary fraction significantly, al- 
though the constraint is weakened slightly because the 
sample size is reduced for higher cutoffs. If the signal-to- 
noise cutoff is reduced to 1, however, the most-probable 
inferred binary fraction jumps from 0.14 to 0.24. This 
is due to the inclusion of a single star, labeled Car-1543, 
with a velocity change At; = 15.7 km/s and Av/a2e — 
3.6. If this star is indeed a binary, then the inferred bi- 
nary fraction is less than 0.42 at 68% CL, and less than 
0.73 at 95% CL; this significantly increases the proba- 
bility that Carina's binary population is consistent with 
Milky Way field binaries, although the inferred binary 
fraction is still low. However, the radial velocity mea- 
surements for this star, 206.6 km/s and 222.3 km/s (com- 
pared to Carina's systemic radial velocity of 223 km/s), 
each have relatively low signal-to-noise, equal to 1.09 and 
1.19 respectively, and thus may be suspect. 

A further complication arises in that, although hori- 
zontal branch stars cannot be included in our analysis, 
there are seven stars in the Carina sample at magnitudes 
My > 20.39 (placing them on or near the horizontal 
branch) that show significant velocity variations, greater 
than 15 km/s and as high as 40 km/s. None of these stars 
have colors that place them in the RR Lyrae instability 
strip, so they are unlikely to be variable stars. Instead, 
four of them lie in the red clump horizontal branch region 
of the color-magnitude diagram, two are clear Milky Way 
stars due to their low velocities, and only one member 
star lies on the red giant branch. However, three of these 
stars have signal-to- noise less than 1.1, and all of them 
have measured velocities with low Tonry-Davis i?-values, 
between 4 and 5 (measurements with i?-values less than 
4 were discarded from the published sample to avoid non- 
Gaussian error due to false-peak selection) . Since the ve- 
locity variations are quite high compared to what might 
be expected from binaries, we find it likely that these 
stars show velocity variations because of non-Gaussian 
error related to their relatively faint magnitudes. While 
the other galaxies in the sample contain relatively few 
stars showing variability at these magnitudes, we note 
that the Carina sample contains 60 multi-epoch stars 
with Aly > 20.3, while the other three galaxies combined 
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binary fraction (B) 

Fig. 6. — Posterior probability distributions in binary fraction 
for tlie Carina sample, assuming different signal-to-noise cuts. The 
results are essentially tiie same for S/N cuts greater than 1.2, but 
the most probable inferred binary fraction increases to ^ 0.24 if 
stars with S/N > 1 are included (solid curve). This is due to the 
inclusion of a single star with signal-to-noise measurements less 
than 1.2. If this star is a binary, the possibility that Carina's 
low apparent binary fraction is a statistical fluke becomes more 
significant, as binary fractions up to 0.73 are allowed to within 
95% confidence limits. 

only have 19 multi-epoch stars in this magnitude range. 
Thus, we consider it unlikely that the velocity variations 
of stars in Carina beyond My > 20.3 are primarily due 
to binary motion. 

8. PERIOD DISTRIBUTION CONSTRAINTS 

In the previous section we derived binary fraction con- 
straints by assuming the period distribution of each 
galaxy is identical to that of Milky Way field binaries, 
which takes the form of a log-normal distribution with 
mean /uiogp = 2.24 and dispersion aiogp = 2.3. Here we 
shall likewise assume the period distribution to have a 
log-normal form, but will allow /xiog p and criog p to vary. 
Thus the likelihood in each star is given by eq. [51 where 
our set of binary parameters V — {/iiogP, criogp}- As in 
section [721 also use magnesium strength information 
to help determine membership and constrain the member 
fraction F. Thus our four parameters being varied are 
the member fraction F and the set of binary parameters 

{B, /XlogP, CTlogp}- 

As in the previous section, since the binary fraction B 
may have any value between and 1, we choose a uniform 
prior in B over this interval. Our prior in the period dis- 
tribution parameters ^,\ogp and criogP, however, requires 
more careful consideration. The prevailing viewpoint is 
that the observed distribution of field binary stars in the 
solar neighborhood is a superposition of populations from 
a wide variety of star-forming environments with differ- 
ent period distributions; this is supported by the fact 
that several clusters seem to exhibit period distributions 
that may be more peaked than those observed in the field 
(jBrandner fc Koehleill998l:IScallv et"allll9 99l. However, 
the period distributions in these clusters still range over 
multiple decades of period. To be conservative, we shall 
consider a flat prior in the mean log-period fiiogp over 
the range [0,4], where the lower limit in this range cor- 
responds to a mean period of one year; as we shall see 
shortly, a binary population with such a short mean pe- 



riod would produce velocity variations in excess of what 
is observed in each galaxy, unless the binary fraction is 
quite small. We will also choose a conservative prior over 
the period distribution width uiogp over the range [1,3]. 
The lower limit ciogp = 1-0 is narrower than any open 
cluster yet observed, and the field binary populations in 
larger dwarf galaxies might be expected to be broader 
than that of clusters; on the other hand, given the rela- 
tively small size of dwarf galaxies, it would seem unlikely 
for the period distribution to be considerably broader 
than that of the Milky Way. In practice, the range of 
the priors will only be relevant when marginalizing over 
the allowed range of ciog p to obtain probability distribu- 
tions in B and /xiogp, as a more conservative prior would 
produce weaker constraints in these parameters. 

When attempting to constrain the period distribution 
of a binary population, it becomes desirable to include as 
many repeat measurements as possible for each star. Un- 
fortunately however, calculating the binary likelihoods 
becomes computationally expensive for more than two 
epochs, because the binary tail in the likelhood becomes 
spread out over the velocity space {Avi, Av2, ■ • •) and 
thus a large number of Monte Carlo points are required 
to generate a smooth likelihood, especially for stars with 
large velocity changes in which case the likelihood is quite 
small. By running our code in parallel over a large num- 
ber of processors, we can generate likelihoods for stars 
with measurements over as many as four epochs. Only 
the Sextans sample contains stars with measurements at 
five epochs or more, albeit only a small number. For 
these stars, we discard any epochs occurring after the 
fourth measurement. 

For each galaxy, we calculate the likelihoods for each 
star in a table of values over the period distribution 
parameters and subsequently interpolate to find the bi- 
nary likelhood for any combination of binary parameters 
MiogP, and (Tiogp. After generating the likelihoods for 
each star, we obtain joint posterior probability distribu- 
tions in /^logp vs. B, and /iiogp vs. ciogp by marginal- 
izing over the other parameters using a nested sampling 
routine. Contour plots of the resulting posteriors for each 
galaxy are shown in Figure [71 For each galaxy, we find 
that the resulting posteriors are degenerate in the three 
binary parameters — for example, a population with a low 
binary fraction and short mean period, and a popula- 
tion with a high binary fraction and long mean period, 
are equally allowed by the data. These degeneracies re- 
sult from the fact that a large fraction of short-period 
binaries can be obtained by having a large binary frac- 
tion, a short mean period (small /xiogp), or a wide period 
distribution (large criogp). The degeneracy between the 
binary fraction and the period dist ribution parame ters 
has been investigated in detail in Mi nor et al.l (|2010( ) . In 
the next section we shall discuss what would be required 
to break this degeneracy and obtain strong independent 
constraints on the binary fraction and period distribution 
parameters. 

By way of comparison, Milky Way field binaries in the 
solar neighborhood have binary parameters B — 0.5, 
MiogP — 2.24, CTiogp = 2.3. This set of parameters 
lies squarely within the allowed parameter region for the 
Fornax, Sculptor, and Sextans dSphs, indicating that a 
Milky- Way like binary population is compatible with the 
data in those galaxies. On the other hand, the Milky Way 
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(b) Fornax 




(d) Carina 




(f) Sculptor 



(g) Sextans 
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Fig. 7. — Joint posteriors in binary fraction vs. mean log-period (left), and dispersion of log-period vs. mean log-period (right) for each 
galaxy. A flat prior is assumed for each parameter. Milky Way field binaries have a binary fraction 0.5 for solar-type stars, with a mean 
log-period /^logp = 2.24 and log-dispersion of periods o"iogP = 2.3. These values lie squarely within the allowed parameter space for the 
Fornax, Sculptor, and Sextans dSphs, but are outliers for Carina. Compared to Milky Way field binaries. Carina's binary population likely 
has either a smaller binary fraction, a longer mean period, or both. 
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Fig. 8. — Joint posteriors in binary fraction vs. mean log-period l |8(a)[ l, and mean log-period vs. dispersion of log-period l |8(b)[ l for a 
simulated sample of 1500 stars with measurements over 2 epochs, and a measurement error of 1.5 km/s. Even in this large sample, there 
is a clear degeneracy between the binary fraction and mean log-period, and also between the period distribution parameters 



parameter set is an outlier for Carina, falling on the edge 
of its allowed parameter space. Compared to Milky Way 
field binaries, Carina's binary population likely has ei- 
ther a smaller binary fraction (B < 0.5), a longer mean 
period (niogP > 2.3), or a combination of these. 

For each galaxy, we sec that parameter regions where 
the mean period is short and the binary fraction is suf- 
ficiently high are ruled out, because such binary popula- 
tions would produce greater velocity variations than are 
observed in the data. On the other hand, regions where 
the mean period is long and binary fraction is low are 
also ruled out, since they would produce fewer velocity 
varations than are observed. Once again the Carina dSph 
stands out, in that it allows small binary fractions even 
if the mean period is relatively long (up to 10,000 years). 
This reflects the lack of velocity variations observed in 
Carina. However, in the region where the mean period 
is longer than that of the Milky Way (/iiogP > 2.24), the 
binary fraction is essentially unconstrained, for a simple 
reason: even if the binary fraction is high, binary motion 
can be "hidden" if the width of the period distribution is 
small, restricting all binaries to long periods; conversely, 
if the period distribution is broadened, the binary frac- 
tion can be reduced so that very few short-period bi- 
naries would be observed. Thus there is a degeneracy 
between binary fraction and the width of the period dis- 
tribution that prevents us from constraining the binary 
parameters if the mean period is assumed to be similar 
to the Milky Way or longer. We can say, however, that 
if the mean period is of order 30 years (/iiogP ~ 1-5 or 
shorter, the binary fraction must be smaller than 0.4 re- 
gardless of the width of the period distribution, and is 
likely smaller than 0.2. We emphasize again that these 
results are somewhat prior-dependent, in that a more 
conservative prior on crjog p would widen the allowed pa- 
rameter region somewhat. 

9. BINARY CONSTRAINTS IN SIMULATED DATA SETS 

9.1. Degeneracy between binary fraction and period 
distribution parameters 

Since the likelihood method demonstrated in previous 
sections can constrain the parameters characterizing a 
binary population, it is natural to ask: how large a sam- 
ple, and how many epochs, would be required to obtain 



strong and independent constraints on the binary frac- 
tion and period distribution parameters? To give some 
idea of this, we first simulate a binary population with 
a binary fraction B = 0.5, and period distribution pa- 
rameters identical to that of Milky Way field binaries 
(MiogP = 2.24, CTiogp = 2.3). We assume a measurement 
error of 1.5 km/s, which is attainable from multi-object 
spectrographs, especially if measurements are averaged 
together over short timescales. After using a Monte 
Carlo simulation to generate mock stellar line-of-sight 
velocities in a large sample of 1500 stars, we then obtain 
posterior probability distributions using the same proce- 
dure as in section[5] Posteriors in B vs. /iiogp and /iiogp 
vs. (Tiogp are p lotted in Figure |S1 



Figure |8(a) shows that while the allowed parameter 
space in B and fiiogp has narrowed considerably com- 
pared to the data (Figure [7]) , there is a clear degeneracy 
between these parameters: a high binary fraction and 
long mean period can produce the same observed binary 
variation, to within the measurement errors, as a low 
binary fraction and short mean period. A similar de- 
generacy exists between the binary fraction and spread 
of periods criogp. The reason for th is degeneracy, which 
is explored in detail in IMinor et al.l (2010), is that only 
binaries with periods in the short-period tail of the pe- 
riod distribution have observable variations beyond the 
measurement error. Thus the shape of the period dis- 
tribution is not well constrained, as is evident in Figure 



8(b)| and adjusting the mean period or spread of periods 
can be largely compensated for by changing the binary 
fraction. We find that this degeneracy is still present 
even if measurements are taken at 3 or 4 epochs, indi- 
cating that the measurement error is simply too large 
to constrain the shape of the period distribution inde- 
pendently of binary fraction. We note, however, that if 
the mean period of a binary population is shorter than 
that of Milky Way binaries, then the period distribution 
can be constrained more easily since velocity variation 
is observable for binaries with periods covering a larger 
portion of the period distribution. 

9.2. Breaking the degeneracy between binary fraction 
and period distribution parameters in a simulated 
dataset 
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Fig. 9. — Marginal posterior probability distributions in each binary parameter, for a simulated star population with binary fraction 
B = 0.5, mean period of 10 years {fiiogp = 1), and log-spread of periods ui^^p = 1.5. Each star has measurements taken at 4 epochs, each 
spaced 1 year apart, and the measurement error is 0.5 km/s for each measurement. As the sample size is increased, the degeracy between 
the binary parameters becomes weaker and the posteriors peak around the true values in each parameter. 



Before we tackle the question of what would be re- 
quired to constrain the binary fraction and period distri- 
bution parameters independently in dSphs, it would be 
reassuring to verify that this is even possible in principle 
using our likelihood approach. To show this, we simu- 
late a binary population with a mean period of 10 years 
(WogP = 1), a spread of periods aiogp = 1.5, and a bi- 
nary fraction B = 0.5. We choose a measurement error 
of 0.5 km/s and generate simulated samples of 500, 1000, 
2000, and 5000 stars with this binary population. Since 
binary periods are better constrained with a large num- 
ber of epochs, we assume measurements taken at four 
epochs (the largest number allowable at present; com- 
putations become prohibitively expensive beyond four 
epochs), each spaced a year apart. 

From the simulated line-of-sight velocity data, we ob- 
tain posterior probability distributions via a nested sam- 
pling routine in each parameter B, /iiogp, criogP, as 
shown in FigurelHl Even for a sample of 500 stars, the de- 
generacy between binary fraction B and mean log-period 
MiogP is being broken, as evidenced by the peaks ap- 
pearing in the posteriors near their correct values. The 
width of the period distribution uiogp is more difficult 
to constrain, and only becomes well-constrained beyond 
1000-2000 stars. The degeneracy between all three pa- 
rameters is clearly broken beyond 2000 stars. Note that 
once the degeneracy is broken, the error in the maximum- 
probability values in each parameter (i.e. the difference 
between the peak value and the "true" value) is con- 
siderably smaller than the width of the posterior distri- 
butions, which suggests that the Bayesian error may be 
overly conservative in this case. 

How important is the number of epochs for constrain- 
ing the period distribution parameters? In Figure [TU] we 
plot posterior distributions in B for 2, 3, and 4-epoch 
samples of 2000 and 5000 stars. If measurements are 
taken at only two epochs, the degeneracy is not broken 
even for a 5000-star sample. Likewise, if measurements 
are taken at three epochs, the degeneracy is not broken 
for a 2000-star sample, although it is broken for a sample 
as large as 5000 stars, as evidenced by the peak appearing 
in the posterior. However, the degeneracy is much easier 
to break in a four-epoch sample; the binary fraction is 
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Fig. 10. — Marginal posterior probability distributions in binary 
fraction, for the same simulated star population as in Figure |9] 
with a binary fraction B = 0.5. We consider samples of 2000 stars 
with 3- or 4-epoch measurements (solid curves) and 5000 stars 
with 2-, 3-, and 4-epoch measurements (dashed, dot-dashed, and 
dotted curves respectively). Note that the binary fraction cannot 
be recovered with only 2 epochs even in a 5000-star sample, since 
it is still degenerate with the period distribution. In a 2000-star 
sample, four epochs are required to constrain the binary fraction 
independently of the period distribution. 

better constrained in a four-epoch sample of 2000 stars 
compared to a three-epoch sample of 5000 stars. We thus 
conclude that if the measurement error is small enough 
to allow the period distribution to be constrained, 3- or 
4-epoch samples of at least several hundred (and pos- 
sibly thousands) of stars are required to obtain strong 
and independent constraints on the binary parameters — 
in other words, to break the degeneracy between the bi- 
nary fraction and period distribution parameters. 

9.3. Required measurement error for obtaining 
independent constraints on binary parameters 

For a binary population with a given mean period, 
what measurement error would be required to constrain 
the period distribution independently of the binary frac- 
tion? In principle, even if only binaries whose periods 
lie on the short-period tail of the period distribution 
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have observable velocity changes, the shape of the pe- 
riod distribution could be inferred if a sufficient number 
of binaries are observed whose periods are well-known. 
This only holds to the extent that the assumption of a 
log-normal period distribution is accurate, so that the 
shape of the period distribution can be extrapolated to 
longer periods. Unfortunately however, with only 2-4 
epochs, binary periods of individual stars cannot be de- 
termined with great accuracy. Furthermore, the likeli- 
hood approach is not computationally feasible at present 
with more than 4 epochs, because of the enormous num- 
ber of Monte Carlo points required to generate smooth 
likelihoods. It is thus reasonable to expect that to con- 
strain the period distribution well using the likelihood 
approach, observing only binaries on the short-period 
tail of the period distribution is not sufficient. We make 
the assumption that binary periods nearly as long as the 
mean period (~ 180 years for Milky Way field binaries) 
should be observable, i.e. should have observable velocity 
variations, in order to constrain the period distribution 
parameters well. We can get an approximate idea of the 
measurement error required for this by considering typ- 
ical velocity amplitudes associated with b inaries with a 
given orbital period. In lMinor et al.l (j2010[ ). a formula is 
derived for the velocity amplitude of a binary as a func- 
tion of period, which is averaged over the possible orien- 
tations, mass ratios, and eccentricities of binary systems 
(assuming Milky- Way like PDFs in these parameters): 



(17) 



where M is the mass of the primary star. If we assume 
M w 0.8Mq which is typical of old stellar populations 
in dwarf spheroidals, we find that binaries with periods 
longer than a few decades will yield velocity amplitudes 
less than 2 km/s, and thus would be difficult to distin- 
guish from a measurement error of 1-2 km/s. Binaries 
with a period of 180 years would have velocity ampli- 
tudes of « 0.9 km/s on average. However, there is a 
further difficulty: if measurements are taken over the 
course of a few years, binaries of such long periods would 
traverse only a small fraction of their orbits and thus 
would exhibit velocity changes much smaller than their 
overall velocity amplitude. As a rough approximation, if 
we assume a circular orbit, it can be shown that the rms 
change in line-of-sight velocity after a time At is given 
by the simple equation 
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where we have averaged over all possible orbital angles. 
In an approximate sense, velocity changes will be ob- 
servable if they are larger than the measurement error 
between two pairs of velocity measurements. If the mea- 
surement error am is roughly the same for either mea- 
surement, then the error in the velocity difference Av is 
o'2e = V^o-jn- Usiug this together with eqs. [T71 and [T51 
we find the measurement error must be smaller than 
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Fig. 11. — Approximate measurement error required to observe 
velocity variation in binaries of a given orbital period, over a 
timescale of 1 year (solid curve) and 5 years (dotted curve) . As can 
be seen from the dotted curve, to observe orbital periods compara- 
ble to the mean period of Milky Way field binaries (logjQ P = 2.24), 
a measurement error on the order of 0.1 km/s would be required 
over a timescale of several years. 

This is the maximum (approximate) measurement er- 
ror with which one can observe binary motion for a bi- 
nary with period P over a timescale At. In Figure [TT] 
we plot the required measurement error as a function of 
log of period over timescales of 1 year and 5 years. Note 
that the required measurement error is very small near 
periods for which the timescale At is an integer multiple 
of one-half the period, since the velocity change is typi- 
cally small for these timescales. However if the period is 
long enough, the required measurement error decreases 
monotonically. For binaries with periods equal to the 
Milky Way's mean period (logP = 2.24), however, the 
maximum measurement error for observing binary mo- 
tion over 5 years is ~0.1 km/s, and much smaller over 
a timescale of 1 year. We can therefore expect that if 
a binary population has a mean period similar to that 
of Milky Way field binaries, it would be extremely dif- 
ficult to constrain its period distribution parameters in- 
dependently of binary fraction if the measurement error 
is significantly larger than 0.1 km/s, unless the measure- 
ments are taken over timescales considerably longer than 
5 years. The prospects are improved if the mean period 
is shorter; for example if the mean period is 10 years. 
Figure [TT] suggests that the period distribution can be 
constrained independently of binary fraction if the mea- 
surement error is less than roughly 1 km/s. 

It should be emphasized that eq. [19] gives only an ap- 
proximate idea of the measurement error required to 
break the degeneracy between binary fraction and the 
period distribution parameters; it says nothing about 
the number of stars and epochs that would be required 
to break this degeneracy, even if the measurement error 
is small enough. As we explored in section 19.21 for the 
case of a relatively short mean period, at least three or 
four epochs and several hundred up to several thousands 
of stars would be required, but this requirement would 
likely vary depending on the actual period distribution 
and measurement errors in question. 

Finally, we note that the task of constraining the bi- 
nary period distribution would be considerably easier 



17 



if the binary fraction were known a priori. This is 
possible in globular clusters, where the binary fraction 
can be con strained by photometry (jSoUima et alj [20071 
iSoUima et" al. 2012). I3ecause the stars in the cluster have 
nearly identical ages and metallicities, the main sequence 
is well-defined and binaries will lie off the main-sequence 
for single stars due to the combined light of the primary 
and secondary star. After the binary fraction is con- 
strained from photometry, this constraint could be in- 
cluded as a prior in a likelihood analysis following a radial 
velocity variability survey to constrain the period distri- 
bution in the cluster. Unfortunately in dwarf spheroidal 
galaxies, constraining the binary fraction through pho- 
tometry is a much more difficult task due to the compli- 
cated distribution of stellar ages and metallicities in each 
galaxy. 

10. DISCUSSION 

We demonstrated in section \T7]\ that the Carina dSph 
sample is nearly devoid of short-period binaries, with 
a binary fraction less than 0.5 to within 90% confi- 
dence limits if a Milky Way-like period distribution is 
assumed. While the reason for Carina's apparent short- 
age of close binaries is unclear, several factors suggest 
themselves. Photometric stu dies show Carina to have a 
bursty star formation hist ory ()Smecker-Hane et"al]|1996l : 
iHurlev-Keller et al.l 11998 ^ , with two predominant star 
populations: approximately 25% of its stars are older 
than 10 Gyr, while the remaining majority are interme- 
diate age stars with ages in the range of 4-7 Gyr. The 
question is therefore, what initial conditions during its 
relatively short period of star formation might have pre- 
vented close binaries from forming? Potential factors 
are the initial tempe rature of th e star-forming clouds 
prior to collapse fSterzik et al.l f2003). as well as the pres- 
ence of magnetic fields and strong radiative feedback, 
whic h can all play a role in inhibiting binary star forma- 
tion ([Price fc Bate! 1201 Of ). While Carina's low metaUic- 
ity compared to the Milky Way might conceivably affect 
radiative feedback and hence suppress binary star forma- 
tion, we note that the Sculptor and Sextans dSphs also 
have relatively low mean metallicities, and their best-fit 
binary fractions are both higher than 0.5. Thus, its low 
metallicity alone is unlikely to explain the apparent lack 
of binaries in the Carina dSph. 

It is worth noting that a recent spectroscopic study 
of metal abundances in Carina p^emaslc ct al. 2012) 
showed little evidence of iron enrichment from type la 
supernovae during its intermediate period of star forma- 
tion, even though a few Gyr should have been a long 
enough timescale for SNe la to contribute iron to the 
interstellar medium while star formation was still tak- 
ing place. While there is still some uncertainty in the 
timescale for star formation of Carina's intermediate age 
population, this lack of evidence for SNe la would be 
expected if Carina is indeed deficient in short-period bi- 
naries, as the present study sugg ests. 

As we cautioned in section [731 however, the possibility 
that Carina's apparently low fraction of binaries may be 
a statistical coincidence increases somewhat if the star 
Car- 1543 is a binary, although the best-fit binary frac- 
tion is still low. However, this star's observed velocity 
variation may be explained by low signal-to-noise mea- 
surements (S/N ~1.1). Furthermore, stars near the hor- 



izontal branch region also exhibit velocity variations, al- 
though most or all of these are probably explained by 
poor spectrum cross-correlation fits related to their faint 
magnitudes. 

11. CONCLUSIONS 

We have developed a general methodology for con- 
straining properties of binary star populations from 
multi-epoch line-of-sight velocity measurements in dwarf 
galaxies. This method has been applied to data fr om the 
Magellan/MMFS sample of Wal ker et aP ()2009al ) in the 
Carina, Fornax, Sculptor, and Sextans dwarf spheroidal 
galaxies to find constraints in their binary fraction and 
period distribution parameters. To obtain the best possi- 
ble binary constraints, we have also re-derived the mea- 
surement eiTors_Jn_the samp le by extending the error 
model of i Walker et al.l (|2007f) to account for binary or- 
bital motion. The best-fit binary fractions in each galaxy, 
assuming a Milky Way-like period distribution, are listed 
in Table [5J with probability distributions in binary frac- 
tion plotted in Figure HJ more general probability distri- 
butions in the period distribution parameters and binary 
fraction are plotted in Figure [71 We conclude with the 
following points: 

1. If a Milky Way- like period distribution is assumed 
in each galaxy, the Fornax, Sculptor, and Sextans dSph 
galaxies have binary fractions that are roughly consis- 
tent with that of Milky Way field binaries, whereas the 
Carina dSph is apparently deficient in binaries compared 
to the Milky Way field. Carina's inferred binary fraction 
is less than 0.5 at the 90% confidence level, with a best 
fit value of 0.14lo;o5; thus a Milky Way-like binary frac- 
tion of « 0.5 in Carina is statistically unlikely. Relaxing 
the assumption of a Milky Way-like period distribution, 
the lack of observed binary velocity variation in Carina 
could be explained by either a small binary fraction, or a 
long mean period, compared to the other galaxies in the 
sample (see Figure [7]). 

2. With the exception of Carina, the published 
measurement errors in the Magellan/MMFS sample of 
Walk er et al.l (|2009aD underestimate the Gaussian mea- 
surement error somewhat, as can be seen in Table[TJ This 
is most striking in the case of the Fornax dSph, where 
the median measurement error is larger than the pub- 
lished value by a factor of ~ 55%. While this difference 
of « 0.6 km/s in the median error is small compared 
to Fornax's velocity dispersion of ~ 12 km/s, the extra 
measurement error should be taken into consideration 
when applying mass models to Fornax. This may be 
particularly of concern if higher moments in the veloc- 
ity distribution are used, since stars lying on the tail of 
Fornax's velocity distribution may have unaccounted-for 
measurement error if the published errors are used. 

3. While multi-epoch surveys can in principle produce 
independent constraints on a galaxy's binary fraction and 
the period distribution parameters, these parameters are 
unfortunately degenerate given the present measurement 
errors. For example, a binary population with either a 
low binary fraction and short mean period, or a high bi- 
nary fraction and long mean period, can both produce 
the same velocity variations observed in the data equally 
well. If a galaxy's mean binary period is roughly compa- 
rable to that of the Milky Way, strong independent con- 
straints on the binary parameters can only be obtained 
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with a measurement error of order ~0.1 km/s or smaller, 
obtainable only by a high-resolution spectrograph. In 
the near future, the best possible binary constraints may 
be obtained in globular clusters where the binary frac- 
tion can be estimated independently through photome- 
try, followed by a multi-epoch radial velocity survey to 
constrain the binary period distribution. 

4. Although the degeneracy between the binary 
parameters cannot be broken with measurement er- 
rors presently obtainable by multi-object spectrographs, 
meaningful comparisons between different binary popu- 
lations can still be made. This is demonstrated in Fig- 
ure [T] although independent constraints on the binary 
fraction and mean period cannot be obtained, it is evi- 
dent that the Carina dSph occupies a different region of 
parameter space compared to the other galaxies in the 
sample. Larger multi-epoch surveys in dwarf galaxies will 
narrow the parameter space further to allow even more 
robust comparisons between galaxies. By comparing the 



allowed parameter space for different galaxies, we can 
ultimately address the question of whether field binary 
populations are universal in nature, or vary depending on 
the initial conditions present during the galaxy's epoch 
of star formation. Large multi-epoch surveys, in concert 
with statistical methods like that demonstrated here, will 
provide a powerful test of star formation theories in the 
future. 
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